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Abstract 

In this paper we study a singular stochastic differential equation driven 
by an additive fractional Brownian motion with Hurst parameter H > ^. 
Under some assumptions on the drift, we show that there is a unique 
solution, which has moments of all orders. We also apply the techniques of 
Malliavin calculus to prove that the solution has an absolutely continuous 
law at any time t > 0. 



1 Introduction 

The aim of this paper is to study a stochastic differential equation, driven by 
an additive fractional Brownian motion (fBm) with Hurst parameter H > ^, 
assuming that the drift f{t, x) has a singularity at x = of the form x"", where 

The study of this type of equations is partially motivated by the equation 
satisfied by the d-dimensional fractional Bessel process Rt = \Bf\, d > 2 (see 
Guerra and Nualart [7], and Hu and Nualart [8]): 



Rt = Yt + H{d - 1) / 
Jo 



where the process Yt is equal to a divergence integral, Yt = J* X!f=i 
The process Y is not a one-dimensional fractional Brownian motion (see Eisen- 
baum and Tudor [S] and Hu and Nualart for some results in this direction), 
although it shares with the fBm similar properties of scaling and -^--variation. 
Notice that here the initial condition is zero. 
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We are considering the case where the initial condition xq is strictly positive. 
Using arguments based on fractional calculus inspired by the estimates obtained 
by Hu and Nualart in [S] , we will show that there exist a unique global solution 
which has moments of all orders, and even negative moments, in the particular 
case f{t,x) = Kx~^, if t is small enough. We will also show that the solution 
has an absolutely continuous law with respect to the Lebesgue measure, using 
the techniques of Malliavin calculus for the fractional Brownian motion. As an 
application we obtain the existence of a unique solution with moments of all 
orders for a fractional version of the CIR model in mathematical finance ([3]), 
which is a singular stochastic differential equation driven by fractional Brownian 
motion with the diffusion coefficient being ^/x. 

The paper is organized as follows. In the first section wc will consider the case 
of a deterministic differential equation driven by a Holder continuous function, 
and with singular drift. The case of the factional Brownian motion is developed 
in Section 3. 



2 Singular equations driven by rough paths 

For any s < t, C{[s,t]) denotes the Banach space of continuous functions 
equipped with the supremum norm, and C^{[s,t]) denotes the space of Holder 
continuous functions of order /3 on [s,t\. For any x € C([s,t]) we put 

l|a;||s,t,oo = sup{|a;(r)|,s < r < t} , 

and if a; e C^{[s,t]) we put 

,, ,, , \x(u) — x(v)\ -. 

xU,t,/3 = sup{^^ yi^s<u,v<t} . 

\u — v\'^ 

Fix f3 G (1/2, 1). Suppose that (p : M.^ ^ K is a function such that (p{0) = 0, and 
(fi e C^([0,r]) for all T > 0. Consider the following deterministic differential 
equation driven by the rough path (p 

Xt^xo+ / /(s, Xs)ds + tpit) , (2.1) 

where xo > is a constant. We are going to impose the following assumptions 
on the coefficient /: 

(i) / : [0, oo) X (0, oo) — > [0, oo) is a nonnegative, continuous function which has 

a continuous partial derivative with respect to x such that dxf{t, x) < Q 
for aU t > 0, X > 0. 

(ii) There exists xi > and a > ^ — 1 such that f{t,x) > g{t)x~", for all 

t > and X € (0, xi), where g{t) is a nonnegative continuous function 
with g{t) > for aU t > 0. 
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(iii) f{t, x) < h(t) (l + i) for all t > and a; > 0, where h{t) is a certain 
nonnegative locally bounded function. 

Theorem 2.1 Under the assumptions (i)-(ii), there exists a unique solution Xt 
to equation 112. 1\) such that Xt > on [0,oo). 

Proof It is easy to see that there exists a continuous local solution xt to 
equation (|2.ip on some interval [0, T), where T satisfies T = inf{t > : xt = 0}. 
Then it suffices to show that T = oo. Suppose that T < oo. Then, then 
a:* — > 0, as 1 1 T. Since ip e C^{[0,T]), there exists a constant C > 0, such that 
\<f{t) — (p{s)\ < C\t — s\^, for all s,t e [0, T]. Since Xt satisfies the equation 
dllT]), for aU t e [0, T] we have 

= XT = Xt+J^ f{s,Xs)ds + ip{T)-ip{t). 
Since f{s,Xs) is positive, for all t G [0,T] we have 



a;t < a;t + / /(s, = (^(t) - (^(T) < C(T - t) 



From the assumption (ii), there exist to G (0,T) and a constant iiT > 0, such 
that g{t) > K and Xs € {0,xi) for aU t e [to,T). Then, for all t e [io,T) we 
have 

X git) K K 

f{t,xt) > — > > ^ . 

^ ^ - Xt" - Xt" - C"{T^ tf^ 

Consequently, for all i G [to , T) we obtain 

K{T-tf ^13 _ ^ r f{s,xs)ds <C{T~tf, 



C«(l-a/3) 7t C"(T-s)"'3 

which is a contradiction because 1 — aP — (3 < and t can be arbitrarily close 
to T. Therefore, T = oo. This proves the existence of the solution for all t. 

Now we show the uniqueness. If xi^t and X2,t are two positive solutions to 
equation (|2.ip . then 

Xl,t-X2,t^ / [f{s,Xi^s) - fis,X2,s)]ds . 

Jo 

Because dxf{t, a;) < for alH > 0, x > 0, we deduce 

{xi,t - X2,tY = 2 / {xi^s - a;2,s)[/(s,a;i,s) - f{s,X2,8)]ds < 0. 
Jo 

So Xi^t = X2,t- 

Thus we conclude that there exists a unique solution xt to the equation (j2.ip 
such that Xt > on [0, cx)). ■ 
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Remark 2.2 From the continuity of Xt and f{t,x) and the Holder continuity 
of (p{t), we obtain that for any T > 0, x £ C'^([0, T]). 

The next result provides an estimate on the supremum norm of the solution 
in terms of the Holder norm of the driving function ip. 

Theorem 2.3 Let the assumptions (i)-(iii) be satisfied. If xt is the solution to 
equation h2.1\) . then for any 7 > 2, and for any T > 0, 

l|a^llo,T,oo<^?l,7,/3,T(|xo| + l)cxp|c2,^,/3,T(^l + ||<^|lo^)| , (2-2) 

where C'l^-^^p^T and C2,~i.fj,T are constants depending on ll^^lloToo '^'^^ 

Proof Fix a time interval [0,T]. Let yt — x] . Then the chain rule applied to 
x] yields 

/■* i i-i /■* i-i 

yt = x'^+l f{s,y2)ys ''ds + j y^ ''dip{s). 



(2.3) 

The second integral in (|2.3p is a Riemann-Stieltjes integral (see Young [H]). 
From Assumption (iii), we have 



\yt-ys\ = 7 



- 1 — - / 1-- 

f{u,yZ)yu ''du+ yu '' dif{u) 



1- 

yu 



yu 



du + ^ 



1-1 



yu '' d(p{u] 



where Kt — supjgjg t] h(t). Since 7 > 2, we have 

nt 



1-- 



yu 



du < 



II2/IIL.S0 + \\y\\l.t,lc 



{t-s) 



(2.4) 



(2.5) 



Since a > — 1, we have a > 0.(3 > \ — (3. Thus 1 — a < /?. From Remark 1.1, 
we know that y G C^([0, T]), for any T > 0. A fractional integration by parts 
formula (see Zahle [17]) yields 



yl - dip{u) = {-ly- Df^u ^Dl-"ipt_{u)du, 



(2.6) 



where (pt^{u) = ip{u) — fit), and _D"_^ and D]_°' denote the left and right-sided 
fractional derivatives of orders a and 1 — a, respectively (see [13); defined by 



D^+yu 



1-- 
yu 



r(l - a) I (u- s)" 



1-- 
yu 



{u - r)°+i 



-dr 



(2.7) 



and 



t}^(^P^Hl-o.)f^P^dr\. (2.8) 
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From (|2.7p . and using the Holder continuity of y we obtain 



< C {\\y\\\rX{u - sy- + \\yt:^Au - , (2.9) 



where and in what follows, C denotes a generic constant depending on a, (3 and 
T. On the other hand, from (|2.8p we have 

(2.10) 



IA'r>t-(zi)| < q|^||o,T,/5(t-<+''-^ 

Substituting ([23]) and (|2l0l) into ([SH) yields 



< C (bll',;t(« - =s)-" + \\y\\l~^p{u - .)^(i-^)-") 



x||^|lox,3(i-u)"+''-'dw 

< C'||<^||o,T,,3 



X (||y|li;i(t - sf + ||y||,';|(t - s)^(2"^)) . (2.11) 
Substituting ((2lI1) and ((23)) into ((2^ we obtain 



ll2/lls,t.So + llylls.t,; 



(t-s) + C7|l(^|lo,r,/3 



X (ll2/||'i(i-s)^ + ll2/||',;J(t-5)''''-^^) 



Consequently, using the estimate x i < 1 + x for all a; > 0, we obtain 



\\y\\s,t,p < Ktj 



1_2 i_i 

\\y\\s,t.lo + \\y\\s,t'^ 



X {\\yCi + {i + \\y\UtMt-sr^'~^^ 

which implies 



liy|lL;oo + llyl'' 



s,t, oo 
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Suppose that A satisfies 



A<U7^ ■ (2-12) 



.2C7||(p||o,T,/3 

Then for all s,t E [0, T], s < t, such that < — s < A, we have 



i_2 i_j 

' \\y\\ - 



I s,t. oc ~ 1 1 y 1 1 s,t. oc 

and this implies 



i-i 1- 

' Ms, 



I s.t. oc 



(t-s) 



< \ys\+2KT-/ 

+2CjMo,TAy\QUi-^'^)^ + 

Using again the inequality < 1 + x for all x > and a G (0, 1), we have 

\\y\\s,t,oc= < \ys\+4:KT-f{l + \\y\\s,t.,oo){t-s) 

+2C7||(^||o^T./3 (1 + IIz/IIm.oo) {t - + {t- sf, 

which can be written as 

\\y\\s.t,oo (1 - 2CjMQ,TAt - s)'' - 4ifT7(t - s)) 
< \ys\+AKTj{t-s) + 2it-sf. (2.13) 

Now we choose A such that 

^"(2C^Mk^) ^ (16XT7) ^ {scjMIo.t.p) ■ ^^■^^^ 

Then, for all s, i e [0, T] , s < such that < — s < A, the inequality (|2.13p implies 

||y|lM,oo<2|y,|+C^,^,T, (2.15) 

where C^,/3,t = SKtjT + AT'^ . Take n = [^] + 1 (where [a] denotes the largest 
integer bounded by a). Divide the interval [0,T] into n subintervals. Applying 
the inequality (|2.15p for s = and i = A, we have for all i g [0, A] 

||2/||o,t,oo <2|yol + C^,/3,T . (2.16) 

Applying the inequality (|2.16p on the intervals [A, 2A], . . . , [(n — 2)A, (n — 
1)A], [{n — 1)A,T] recursively, we obtain 

\\y\\o,T,oo < 2"|yo|+2"-iC^,/3,T + --- + C^,;3,T 

< 2i^]+\\yo\ + C^,0,T) 

< 2T(2C-,\\v\\o,T.0)T^^^ y{16KTl)V{8C'y\\v\\o.T.i,)'^ +1 (^^y^^ + C-y t) ■ 
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Therefore, we obtain 

llxllo.T.oo < Ci,^,f3,T{\M + 1) exp |C^2,7,/3.T [l + | , 

which concludes the proof of the theorem. ■ 

3 Singular equations driven by fBm 

Let {B^,t > 0) be a fractional Brownian motion with Hm-st parameter H G 
(1/2, 1), defined in a complete probability space (fi, JF, P). Namely, (i?f^, t > 0) 
is a mean zero Gaussian process with covariance 

E(i3f ) = Rh{s, t) = ^ (<2« + s'" - \t - s\'") . (3.1) 

We are interested in the following singular stochastic diff'erential equation 

Xt^xo+ f /(s, Xs)ds + Bf , (3.2) 
Jo 

where Xq > 0, and the function f{s,x) has a singularity at a: = and satisfies 
the assumptions (i) to (iii). As an immediate consequence of Theorem 12.31 we 
have the following result. 

Theorem 3.1 Under the assumptions (i)-(iii), there is a unique pathwise so- 
lution X = {Xt,t > 0) to Equation i3.2\) . such that Xt > t almost surely on 
[0, oo) and for any T > 0, ||X||o,t,oo e LP{^), for all p> 0. 

Proof Fix P e i^,H) and T > 0. Applying Theorem [231 we obtain that 
there is a unique pathwise solution X ~ {Xt, t > 0) to Equation p.2p . such that 
Xt > almost surely on [0, oo) and 



|^||oxoo<Ci,^,0,T(|a;o| + l)exp C2,w l + ||B«||o^5^ \ . (3.3) 



If we choose 7 such that 7 > 2^-1 ; then < 2, and by Fernique's theorem 

(see [6]), we obtain 

E(eC||s"llo^)<oo, (3.4) 

for all C > 0, which implies that E(||X||pj,^) < 00 for all p > 1. ■ 

Theorem 13.11 implies the existence of a unique solution to the following 
stochastic differential equation with non Lipschitz diffusion coefficient: 

Yt = yQ+ f f{s,Y,)ds+ [ y/Y,dB^, (3.5) 
Jo Jo 

where tjq is a positive constant and / is a nonnegative continuous function 
satisfying the following conditions: 
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(a) There exists a:i > such that f{t,x) > g{t) for alH > and x € {0,xi), 

where 17 is a continuous function such that g(t) > if t > 0. 

(b) f{t, x) > xdxf{t, x) for all t > and x > 0. 

(c) f{t,x) < h{t){x + 1) for aU t > and a; > 0, where ft, is a nonnegative 

locahy bounded function. 

The term -y/Y^ appears in a fractional version of the CIR process in financial 
mathematics (see |3] ) and cannot be treated directly by the approaches in Lyons 
PU] , Nualart and Ra§canu 12J, since function g{x) = ^Jx does no satisfy the 
usual Lipschitz conditions commonly imposed. We make the change of variables 
Xt = Then, from the chain rule for the Young integral, it follows that 

a positive stochastic process Y — {Yt,t > 0) satisfies p.Sp if and only if Xt 
satisfies the following equation: 

Xt^2^+ f ^i^^ds + . (3.6) 

Let fi{t,x) — 2f{t,x)x^^. Then fi{t,x) satisfies all assumptions (i)-(iii), and 
hence from Theorem 13.11 we know that there exists a unique positive solution 
Xt to equation (13. 6p with all positive moments. So Yt = Xi /A is the unique 
positive solution to Equation ()3.5p . and it has finite moments of all orders. 

The next result states the scaling property of the solution to Equation (|3.2|) , 
when the coefficient /(s, x) satisfies some homogeneity condition on the variable 

X. 

Proposition 3.2 (Scaling Property) We denote by Eq{xo, /) the equation Ii3.2\) . 
Suppose that xo > 0, and f{t,x) satisfies assumptions (i)-(iii), and f{t,x) is 
homogeneous, that is, f{st, yx) = s™y"f{t, x) for some constants m, n. Then, 

the process (^a^ X± ,t > 0^ has the same law as the solution to the Equation 

tjq{a xa,a j). 

Proof For each a > 0, we know that {a~^ B^^, t > 0} is a fractional Brownian 
motion. We denote Xa.t the solution to the following equation: 

Xa,t^x„+ f f{s,Xa,s)ds + a-"Bg. 
Jo 

So {Xt,t > 0) (the solution to Eq{xo, /)) has the same distribution as {Xa^t, i > 0). 
Then 

a"X^i. = a"xo+ [\"f{s,Xa.s)ds+ B^ 
Jo 

t 

a"-^-"'-"" f(r,a" Xa.,^)dr + Bf, 
which implies the result. ■ 



a^xo + 
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As an example, we can consider the function f{t,x) — s^x where a > 
— 1, and 7 > 0. Then, if {Xt, i > 0) is the solution to Equation 

Jo 

(|3.2p . then (^a^X±,t > 0^ has the same law as the solution to the Equation 

Xt = a"xo + a^-"-f^-T-i f s''X-°'ds + Bf . 

Jo 

3.1 Absolute continuity of the law of the solution 

In this subsection we will apply the Malliavin calculus to the solution to Equa- 
tion (|3.2p in order to study the absolute continuity of the law of the solution 
at a fixed time t > 0. We will first make some preliminaries on the Malliavin 
calculus for the fractional Brownian motion, and we refer to Decreusefond and 
Ustiinel [Ij, Nualart [TT] and Saussereau and Nualart [H] for a more complete 
treatment of this topic. 

Fix a time interval [0,T]. Denote by £ the set of real valued step functions 
on [0, T] and let Ti be the Hilbert space defined as the closure of £ with respect 
to the scalar product {l[Q^t],'^[o,s])n = Rnit^s), where Rh is the covariance 
function of the fBm, given in (|3.ip . We know that 



RH{t,s) ^ an 

Jo Jq 



t rs 

2H^ 



\r — u\ dudr 



tAs 

KH{t,r)KH{s,r)dr, 



where Knit, s) ^ chs^^ " jl{u-s)" 2u" idul{s<t} with eg = gf^^-ln,! !^) 
and B denotes the Beta function, and aH = H{2H — 1). In general, for any 
tfijip G £ we have 



T pT 

2H-2, 



{ip,i))n ^ oiH I I \r-u\ ifrTpududr . 
Jo Jo 

The mapping l[o.t]' — ^ can be extended to an isometry between TC and the 
Gaussian space Hi spanned by B^ . We denote this isometry by ip i — > B^{(p). 
We consider the operator K'^ : £ L^(0, T) defined by 

{K*hp){s)^ r ^{t)^{t^s)dt. (3.7) 



Notice that (Xlf (l[o^t]))(s) — iirH(i, s)l[o^t] (s). For any (f^ij} € £ we have 

= {K*hP. K*h^)l'^o.t) - E(i?^(^)B«((/))), (3.8) 
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and provides an isometry between the Hilbert space H and a closed subspace 
of L2([o,r]). We denote Kh : i^([0,r]) ^ Hh := iCff (^^([0, T])) the operator 
defined by {KHh){t) := Knit, s)h{s)ds. The space Hh is the fractional 
version of the Cameron-Martin space. Finally, we denote by Rh = Kh ° : 
Ti, Hh the operator Rnf = Jq Kh{-, s){Kfjip){s)ds. We remark that for any 
(p (z H, Rnf is Holder continuous of order H. 

If we assume that f2 is the canonical probability space Co([0,T]), equipped 
with the Borel cr-field and the probability P is the law of the fBm. Then, the 
injection Rh ■ H embeds H densely into and {^,H, P) is an abstract 
Wiener space in the sense of Gross. In the sequel we will make this assumption 
on the underlying probability space. 

Let S be the space of smooth and cylindrical random variables of the form 

F = /(S«(^i),...,S«(^„)), (3.9) 

where / G Cg^(M") (/ and all its partial derivatives are bounded). For a random 
variable F of the form p.9p we define its Malliavin derivative as the 7i-valued 
random variable 

i—l * 

We denote by D^'^ the Sobolev space defined as the completion of the class S, 
with respect to the norm 



l^lll.2 



E(F^)+E(||Z?F||^; 



The basic criterion for the existence of densities (see Bouleau and Hirsch [T]), 
says that if F G 13^''^, and > almost surely, then the law of F has 

a density with respect to the Lebesgue measure on the real line. Using this 
criterion we can show the following result. 

Theorem 3.3 Suppose that f satisfies the assumptions (i)-(iii). Let Xt be the 
solution to Equation ^3.2fl . Then for any t > 0, Xt £ D^'^. Furthermore, for 
any t > Q the law of Xt is absolutely continuous with respect to the Lebesgue 
measure on R. 

Proof Fix a time interval [0,T], and let /? G {^,H). We want to compute 
the directional derivative {DXt, (p)^^, for some ip E H. The function h — Rnf 
belongs to C''([0,T]) and /iq = 0. Taking into account the embedding given by 
Rh '■ H — >ri mentioned before, we will have 



{DXt,v)^ = ^U=o, (3.10) 



where Xt is the solution to the following equation 

X^^xo+f f{s,Xt)ds + B^ + eht, (3.11) 
Jo 
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te [0,T], where e G [0,1]. 

From Theorem 13. 1[ it follows that there is a constant C independent of e 
such that 

E ( sup |X"|n < C < oo, 

\0<t<T J 

for all p > 1. From equations (|3.2p and (|3.1ip . we deduce 

X^-Xt= f {f{s,Xt)~f{s,X,))ds + eht. (3.12) 
Jo 

By using Taylor expansion, the equation (|3.12p becomes: 

Xt-Xt= f Qs{Xl-Xs)ds + eht, (3.13) 
Jo 

where 9s — dxf{s,Xs + ds{X^ — Xg)) for some 61 between and 1. By using 
p.7p the solution to equation (|3.13p is given by 

X^-Xt = exp (^J^ Qrdr^ d(i?H^)(s) 

= e exp I^J' Ordr^ dKH^{s,u) ^^^^^ ^^^^^^^ 



lo {In """^ {Is ^^"^0 ^^It^^^) iK*Hv){u)du. 



Using p.7p and (|3.8p we can write 



Xt-Xt = {K*H^){u) (^K*H (^exp (^^ O.dr^^ j {u)du 



e{(p, exp ( / Qrdr 

' H 

ft ft / ft 



= eaH J J ip{s) exp Qrdr^ \s — u]'^^ ^duds 

Since dxf{t, x) is continuous and dxf{t, 2:) < for alH > and x > 0, we have 
exp ^ Qrdrj < 1. As a consequence, 

lim ~ = ^ ^(s)exp dxf{r,Xr)dr^\s-u\'^"-^duds 

= (^'P,exp(^J dxf{r,Xr)dr^l[oj]'^ , 

where the limit holds almost surely and in L'^{fl). Then, taking into account 
(Prrg)) . by the resuhs of Sugita [l5], we have Xt € B^'^, and 

DXt - exp (^J dx.f{r,Xr)dr^ l[o,t]. (3.14) 
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Finally, 

ft rt / rt 



X exp (^J d^f{r, Xr)dr^ \s - u\'^"-'^duds > 0. 

■ 

In the next proposition we will show the existence of negative moments for 
the function f{t, x) — where i^T > is a constant. The proof is based again 
on the techniques of Malliavin calculus. 

Proposition 3.4 Let {Xt,t > 0) be the solution to the equation 

Xt = xo+ f ^ds + Bf. (3.15) 

Then, for all p > 1 and t < ^ (p-yi)H ^ have E{X^^) < Xq^ . 

Proof Obviously the function f{t,x) — ^ satisfies all the conditions (i)-(iii). 
From Equation (|3.14p . we have 



So, 



DsXt = - f ^DsXrdr + l[o.t](s). 

Jo 

/■* K 

I?sXf=exp{-/ — dr}l[o,t](s). 

J s 



For any fixed p > 1, we construct the family of functions fftix) — (^^^^y , 
X > 0. Then Lp^ | a;~^, as e | 0. For each e > 0, (/^e is a bounded continuously 
differentiable function and all its derivatives are bounded. 
By the chain rule we obtain, 

V,{Xt) = ^eixo)+ f cp[{Xs) — ds+f ip'^iX,)dB^ 

Jo Xs Jq 

1 



'P / 7 TF^r-rdB" (3.16) 



Then, Proposition 5.3.2 in \\A^ implies that 
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where 5 is the divergence operator with respect to fractional Brownian motion. 
Substituting ((3l7)) into ((37l6)) yields 

Fix some t, such that K - {p+l)Ht^"^^ > 0, that is t < (^^q^jp-) . Taking 
expectation on above inequality, we have 

Let e tends to 0. Applying monotone convergence theorem, for any fixed p > 1, 
we obtain 
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